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I. INTRODUCTION 



The experimental data for the semileptonic decays of D mesons are unfor- 
tunately not good enough to clearly determine the dependence of the form 
factors. What is known experimentally, apart from the branching ratios, are 
the form factors at one kinematical point, assuming a pole - type behavior 
for all the form factors. This assumption seems reasonable, but within heavy 
quark effective theory (HQET) the kinematic constrait on the form factors 
at = cannot be satisfied unless a special relation is imposed between 
the pole masses and residues. Recently we have developed a model for the 
semileptonic decays D Vlvi and D — > Plvi, where P and V are light 

= 0~ and 1~ mesons, respectively |T||. This model combines the heavy 
quark effective theory and the chiral Lagrangians. HQET is valid at a small 
recoil momentum [0, |^ and can give definite predictions for heavy to light 
{D ^ V or D ^ P) semileptonic decays in the kinematic region with large 
momentum transfer to the lepton pair. Unfortunately, it cannot predict 
the dependence of the form factors 0, 0]. For these reasons, we have modi- 
fied the Lagrangian for heavy and light pseudoscalar and vector mesons given 
by the HQET and chiral symmetry . Our model |T]] gives a natural expla- 
nation of the pole-type form factors in the whole range, and it determines 
which form factors have a pole - type or a constant behaviour, confirming the 
results of the QCD sum rules analysis [Q. To demonstrate that this model 
works well, we have calculated the decay widths in all measured charm me- 
son semileptonic decays The model parameters were determined by the 
experimental values of two measured semileptonic decay widths. The predic- 
tions of the model are in good agreement with the remaining experimental 
data on semileptonic decays. 

The nonleptonic D meson decays are challenging to understand theoret- 
ically (see e.g. and references therein). The short distance effects are 
now well understood [Q, but the nonperturbative techniques required for 
the evaluation of certain matrix elements are based on approximate mod- 
els. Often the factorization approximation is used (see e.g. p and references 
therein). The amplitude for the nonleptonic weak decay is then considered as 
a sum of the "spectator" contribution and the "annihilation" contribution, 
the direct annihilation of the initial heavy meson. In the determination of 
the "spectator" contribution one uses the knowledge of the hadronic matrix 
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elements calculated in D meson semileptonic decays. Another problem in the 
analysis of nonleptonic D meson decays is the final state interactions (FSI) 
[0, H, ^, [ly, [TT|, |r^. These arise from the interference of different isospin 
states or the presence of intermediate resonances, and both spectator and 
annihilation amplitudes can be affected. The FSI are especially important 
for the annihilation contribution, which can often be successefuUy described 
by the dominance of nearby scalar or pseudoscalar resonances [0, H, ||, |10|] . 
The effective model developed to describe the D V{P)li'i decay widths [|l| 
contains only light vector and pseudoscalar final states and, therefore, is not 
applicable to the annihilation amplitudes. Consequently, in the present paper 
we apply this effective model to analyze only those D PV, D PP, and 
D VV decays in which the annihilation amplitude is absent or negligible. 
Other FSI might arise as a result of elastic or inelastic rescattering. In this 
case, the two body nonleptonic D meson decay amplitudes can be written in 
terms of isospin amplitudes and strong interaction phases |13[ . As usual, we 



assume that the important contributions to FSI are included in these phases. 
In fact, we will avoid the effects of the FSI strong interaction phases by 
considering only the D meson decay modes in which the final state involves 
only a single isospin. Our analysis then includes the decays K*^^^, 
D+ p+K^, D+ A'°7r+, D+ K*°p+, D+ <l>7r+, Df $7r+, 
D+ L)° ^ ^T], D+ p^viv') and u;°r/(V). 

To evaluate the spectator graphs for nonleptonic decays we use the form 
factors for the D ^ V and D ^ P weak decays, calculated for the semilep- 
tonic decays [|l|. This explores how well their particular behavior also 
explains the nonleptonic decay amplitudes. At the same time the analysis 
of the nonleptonic decays enables us to choose between different solutions 
for the model parameters found in the semileptonic decays, determining the 
set of the solutions which are in the best agreement with the experimental 
results for the nonleptonic decay widths. Moreover, we obtain a value for 
the parameter (3, which can not be determined from the semileptonic decay 
alone, but enters in the nonleptonic decays. 

The paper is organized as follows. In Sec. II we present the effective 
Lagrangian for heavy and light pseudoscalar and vector mesons, determined 
by the requirements of HQET and chiral symmetry. In Sec. Ill we present 
the results for the D Vlui, D Plui decays [|1|. In Sec. IV we analyze the 
nonleptonic decay widths. Finally, a short summary of the results is given 
in Sec. V. 
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II. THE HQET AND CHPT LAGRANGIAN FOR D V{P)lu 



We incorporate in our Lagrangian both the heavy flavour SU{2) sym- 
metry, and the SU{3)l x SU{3)r chiral symmetry, spontaneously broken to 
the diagonal SU{3)v |jl|| (and references therein), which can be used for the 
description of heavy and light pseudoscalar and vector mesons. The light 
degrees of freedom are described by the 3x3 Hermitian matrices 



n 



K- 



V2 



TT 



K+ 



V6 



V3 



and 



V2 



V2 

Kf 



(2) 



for the pseudoscalar and vector mesons, respectively. The mass eigenstates 
are defined hj r] = rjs cos 9p — rjosinOp and rj' = rjssinOp + rjo cos 9p, where 
6p = (—20 ± 5)° [jl4[ is the rj — rj' mixing angle. The matrices (|^) and (0) are 
conveniently written in terms of 



.ill, 

u = exp(y) , (3) 



where / is the pseudoscalar decay constant, and 



Pf^ = ^y|PM ' (4) 

where gy = 5.9 is given by the values of the vector masses since we assume 
the exact vector dominance Introducing the vector and axial currents 
'^fj. = ^{u^d^u + ud^u^) and = ^{u'^d^u — ud^u^) and the gauge field 
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tensor F^y{p) = d^j^p^ — d^p^ + [PuiPv] the light meson part of the strong 
Lagrangian can be written as 



+ l^tr[F,^{P)F^''m- (5) 

Both the heavy pseudoscalar and the heavy vector mesons are incorpo- 
rated in the 4x4 matrix 



= 1(1+)^)(d:^^/^-Z}„^5), (6) 

where a = 1, 2, 3 is the SU{'i)v index of the hght flavours and D*^^ and Da an- 
nihilate a spin 1 and spin heavy meson cg^ of velocity respectively. They 
have a mass dimension 3/2 instead of the usual 1, so that the Lagrangian is 
explicitly mass independant in the heavy quark limit rric oo. Defining 

Ha = = (^17'^ + ^175)^(1+)^), (7) 

we can write the leading order strong Lagrangian as 



+ tgTr[Hh7MAnbaHa] + il3Tr[H,v^{V^ - p^baHa] 

+ ^Tr{HbHaHaHb) . (8) 

This Lagrangian contains two unknown parameters, g and f3, which arc not 
determined by symmetry arguments, and must be determined empirically. 
This is the most general even-parity Lagrangian of leading order in the heavy 
quark mass (mg — > oo) and the chiral symmetry limit {rriq and the 
minimal number of derivatives). 

We will also need the odd-parity Lagrangian for the heavy meson sector. 
The lowest order contribution to this Lagrangian is given by 
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LoM = tXTr[H,a^,F^''{p)ahHb] • (9) 

The parameter A is free, but we know that this term is of the order 1/A^ 
with being the chiral perturbation theory scale. 



III. FORM FACTORS IN D ^ V/Puil DECAYS 



For the semileptonic decays the weak Lagrangian is given at the quark 
level by the current - current Fermi interaction 

/:e//(AC = AS = 1) = -^[^7.(1 - l>is'l'{l - 75)c] (10) 

where G^? is the Fermi constant, and s' = scosOc + dsinOc, Oc being the 
Cabibbo angle. 

At the meson level we assume that the weak current transforms as (3^, Ir) 
under chiral SU{2i)l x SU{?))r and is linear in the heavy meson fields. In 
our calculation of the D meson semileptonic decays to leading order in both 
1/M and the chiral expansion we have shown that the weak current is [|[] 



Jl^= \ iaTrYi^{l-^,)H,ul] 

+ a2Tr[7^75i/6t;,(p"-V")6c«L] + -, (H) 

where a = fDy/f^D [i]- The ai term was first considered in p|. We found 
[0] that the a2 gives a contribution of the same order in 1/M and the chiral 
expansion as the term proportional to ai. 

The H ^ V and H ^ P current matrix elements can be quite generally 
written as 



< V(,)(e,p')l(^ - Ar\H{p) >= ^l^e'^-z^eXp'^ 

mH + mv(i) 
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-te .q ^q^A^iq ) + t{mH + mv(i)){e^ -q^)A\'{q ) 



le .q 



mn + mv{i) 



2 2 

niff — m 



and 



Vii) 



Af{q') 



'12) 



mjj — m 



p{i) 



Q,]Fi'\q') 



+ 



m\j — m 



p{i) 



(13) 



where, q = p—p' is the exchanged momentum and the index (i) specifies the 
particular final meson, P or V. In order that these matrix elements be finite 
at = 0, the form factors must satisfy the relations 



2m 



V 



'rriH - my 
2mv 



^2(0) = 



(14) 



(15) 



Fi(0) = Fo(0) . 

and, therefore, are not free parameters. 

In order to extrapolate the amplitude from the zero recoil point to the 
rest of the allowed kinematical region we have made a very simple, physically 
motivated, assumption: the vertices do not change significantly, while the 
propagators of the off-shell heavy mesons are given by the full propagators 
l/(p^ — m^) instead of the HQET propagators l/{2mv ■ /c) With these 
assumptions we are able to incorporate the following features: the HQET 
prediction almost exactly at the maximum g^; a natural explanation for the 
pole-type form factors when appropriate; and predictions of fiat behaviour 
for the form factors Ai and A2, which has been confirmed in the QCD sum- 
rule analysis of 0. 

Finally, we include SU (3) symmetry breaking by using the physical mas- 
ses and decay constants shown in Table 1 of ref. [0]. The decay constants for 
the rj and rj' were taken from |15], for the light vector mesons from and 
for the D mesons from 



r|], |T3 and H 



The relevant form factors for D 
our model illl, are 



V decays defined in ([1^ ) calculated in 
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V{i) 



9v 



K 



V(i) 



rriH'ii) 



,"2 fH'({)l3 



+ ai - 5 —- a2 



m 



V{€) 



m 



H 



mv(i) 



K 



V(i) 



2v/rra^ gv 
-a\- 



i / — 



and 



1 



^2 



72 



(18) 
(19) 



(20) 



where the pole mesons and the constants Kv{i)i which contribute to the 
corresponding processes D PV and D V(i)V(2) are given in 

We determined the three parameters (A, ai, ^2) in [|l| using the three 
measured values of hehcity amplitudes V /Vtot = 0.048 ± 0.004, Tl/Tt = 
1.23 ± 0.13 and r+/r_ = 0.16 ± 0.04 for the process D+ K*°l+ui, taken 
from the Particle Data Group average of all the data |1^. The parameter 



j3 could not be determined from this decay rate, since Ao(?^) cannot be 
observed in the semileptonic decays. 

The model parameters appear linearly in the form factors ([T6|)-(pOD, so 
the polarized decay rates Tq, r+ and r_ are quadratic functions of them. For 
this reason there are 8 sets of solutions for the three parameters (A,ai,a2). It 
was found from the analysis of the strong decays D* Dn and electromag- 
netic decays D* D'-f [|19|, that the parameter A has the same sign as the 
parameter A', which describes the contribution of the magnetic moment of 
the heavy (charm) quark. In the heavy quark limit we have A' = — l/(6mc). 
Assuming that the finite mass effects are not so large as to change the sign, 
we find that A < 0. Therefore only four solutions remain. They are shown 
in Table 1. 

The calculated branching ratios and polarization variables for the other 
semileptonic decays of the type D V are in agreement with all the known 
experimental data |I| (see Table 2). 
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In our approach the form factors for D ^ P decays are given by |]T| 



vTLYi—mpf^-^ I ]^ mH'*{i) / 

where the pole mesons and the constants Kp{i), which contribute to the 
corresponding processes D — > PV and D — > P(i)P(2) are given in in We 
neglected the lepton mass, so the form factor Fq, which multiplies q^, did 
not contribute to the decay width. 

Using the best known experimental branching ratio - B[D^ K^l^ui] = 
(3.68 ± 0.21)% [111], we found two solutions for g: 



SOL. 1 : ^ = ^> = 0.15 ±0.08 , 

SOL. 2 : ^ = = -0.96 ±0.18 . (23) 

The quoted error for is mainly due to the uncertainty in the value fo, while 
the quoted error for g^ is mainly due to the uncertainty in f^* ■ Unfortunately 
we were not able to choose between the two possible solutions for g in (|23|) . 
the branching ratios for D ^ P transitions are presented in Table 3. 



IV. NONLEPTONIC DECAYS 



The effective Hamiltonian for charm decays is given by 

= ^V,,V:^{ai{uT^q,MT''c) + a^iuT ^c){q^T^q,)} (24) 

where Vgg' is an element of the CKM matrix, i and j stand for d ot s quark 
flavours, = 7^(1 — 7^), and ai and 02 are the Wilson coefficients: 

ai = 1.26 ±0.04 aa = -0.51 ±0.05 . (25) 
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These values are taken from |ri|, |I2[ and they are in agreement with the 
next-to-leading order calculation 0. The factorization approach in two body 
nonleptonic decays means one can write the amplitude in the form 

<AB\q,r^qjqkr^c\D> = < A|g,r^g,|0 >< ^Ig^F^clD > 

+ <B\qiT^qj\0><A\qkT^'c\D> 
+ < AB\qir^qj\0><0\qkT^'c\D > . (26) 

In our calculations we take into account only the first two contributions. The 
last one is the annihilation contribution, which is absent or negligible in the 
particular decay modes we consider. In other decays this contribution was 
found to be rather important [|ll|, ^ . It was pointed out in ||, |10 



that the simple dominance by the lightest scalar or pseudoscalar mesons in 
< AB\qir^qj\0 > can not explain the rather large contribution present in 
some of the nonleptonic decays, which we will not consider. Our model [||], 
being rather poor in the number of resonances, is applicable to the analysis 
of the spectator amplitudes, but not the annihilation contributions. 

We will use the following definitions of the light meson and the heavy 
meson couplings < P(p)|j^|0 >= -ifpp^,, < V{p, e*)\jf,\0 >= my /ye*, < 
0\j^\D{P) >= -ifDmnv^, and < 0\j^\D*{e,P) >= imo^fo*^^- Then using 



(|T2D and (|T3D we can write the amplitude for the nonleptonic decay D PV 
processes as 

M{D{p)^PV{e*)) = ^e*-p2mv[-wvKv fpAoiml) 

+ wpKp fv F,{m^y)] (27) 

The factors wy, wp, Ky and Kp are given in 0. 
The D ^ P1P2 decay amplitude is 

G 

M{D{p) P(i)P(2)) = ^[-^^1 fm i^H - ^p{i)) ^o^^("^p(2)) 

- iw2Kp(^2) fp{i) {m^ -ml^2)) Fo^\'mpii))] (28) 

The factors Wi, W2, Kp(i) and Kp{2) are presented in 0. 
Finally, we find the D V(i)V(2) decay amplitude to be 

M{D{p) ^ V(i)(pi,ei),\/(2)(p2,e2)) = (29) 
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+ W2Kv{2) fv{i) mv{i) ei^ 



\^WiKv(l) fv{2) mv(2) ^2, 

+ i{mH + mv{i)) A'i\mv(^2)) ' 





4-Pvi {P + PV2T 



4 ■PV2 {P + PvlT 



(^lu Pa Pip 



The factors wi, W2, Kv(i) and Kv{2) for D V(i)V(2) are given in 0. 

In order to avoid the strong interaction final state effects in the interfer- 
ence between different final isospin states we analyze decays in which the final 
state involves only a single isospin. This occurs when there is an isospin zero 
particle in the final state (a;, r], r^'), or when a final state has the maximal 
third component of the isospin; for example, A'*°7r+, p^K*^, 

D+ K^TT+ and D+ ^ K*'^p+ with |J, J3 >= |3/2,3/2 >). 

Our analysis of semileptonic decays D V{P)lh'i [0] left some ambiguity 
in the choice of the model parameters: there are two values of g, ((?< and g^) 
( |23D and four solutions for the parameters (A, ai, 0:2) (Table 1). The calcu- 
lated nonleptonic decay amplitudes depend on the choice of these parameters. 
However, although the uncertainties of the predictions are quite large, they 
are mostly due to the calculated errors in g^ and g^ (p3|), which is in turn due 
to the uncertainty in f£, and Jd* . The only parameter that is not constrained 
by the semileptonic decay data is the parameter (3 in the form factor Aq, but 
the predictions for the nonleptonic decay rates are not very sensitive to (3. 
From (|27| ) and ([17|) it can easily be seen that (3 appears multiplied by mp in 
the D PV decay width and is only significant for the decays D PV, 
where P is K, t] or t]'. 

First we discuss the results for the decay amplitudes which depend only 
on the form factors Fq and Fi and consequently only on the parameter g; 
namely, D+ K°n+, D+ $7r+, £)+ p^viv'), ^ and 
$7r°. The comparison with the experimental data does not exclude either 
of the values for g, g^ or (?> 0. Next, we summarize the results obtained 
for the decays which depend only on the form factors V, Aq, Ai and A2, 
and consequently only on the parameters (A, ai, 0:2); namely, Dj^ — > $7r"^, 
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D+ <l>p+, ^ •I'pO and D+ K*^p+. The decay D+ <i>n+ depends 
also on the parameter /3, but this dependence is very shght, since the light 
pseudoscalar meson in the final state is a vr. 

The results for all sets are in rather good agreement with the experimental 
data, with the exception of which we do not understand. 

In addition to the above two types of nonleptonic decays, there are two 
measured branching ratios for — »• K*^7r^ and — > p^K^. Their decay 
amplitudes depend on both g and the parameters A, ai, a2. The branching 
ratio for K*°7r+, which is not sensitive to (3 since the vr mass is 

small, excludes the parameter the sets 11 and IV, and prefers g = g> = 
0.15 ±0.08 and the set l._ 

/^From the — ^ p^K^ decay, which has K pseudoscalar meson in the 
final state, one can then estimate the parameter j3. Unfortunately, this decay 
has a considerable experimental error, BR = (6.6 ±2.5)% which results 
in large error in /?: 

^ = 3.5 ±3 . (30) 

The predictions for the branching ratios for the possible decays are pre- 
sented in Table 4 assuming set 1 for A, ai and a2, g = g> = 0.15 ± 0.08 and 
/? = 3.5 ± 3. The quoted errors are due to the uncertainties in the model 
parameters, mainly g. 



VI. SUMMARY 



We have proposed a method to include the light vector meson resonances 
in the weak currents using HQET and CHPT. Instead of the propagators used 
in HQET we have used full propagators for the intermediate heavy meson 
states. In this way we obtain a pole-type behavior of the form factors for the 
matrix elements of the vector currents, and a constant behavior of the form 
factors of the axial current. The calculated branching ratios are in agreement 
with the experimental results. We have predicted the other semileptonic 
decays that have not yet been observed. In addition we have calculated the 
branching ratios for the nonleptonic decay modes D —>■ PV, D P1P2 and 
D ^ V1V2 in which the annihilation contribution is absent or negligible, and 
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in which the final state involves only a single isospin in order to avoid the 
effects of strong interaction phases. Factorization of the matrix elements 
was then assumed and we used the effective model developed to describe 
the semileptonic decays D V{P)lvi to calculate the nonleptonic matrix 
elements. We reproduced the experimental results for branching ratios for 
the K*°-K+, D+ p+K^, D+ $7r+, D+ p+r], D+ K^7r+, 

$p+ and — > K*^p^ decays, albeit within substantial uncertainties. 
We also determined the set of parameters A, ai , 02 and g, which gave the 
best agreement with the experimental results and used this set of parameters 
to estimate the parameter f3 from the branching ratio for p^K^ . We 

then made the predictions for a number of nonleptonic decay rates which 
have not yet been measured. 
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A [GeV-i] 


ai [GeVi/2] 


a2 [GeVi/2] 


Set 1 
Set 2 
Set 3 
Set 4 


-0.34 ±0.07 
-0.34 ±0.07 
-0.74 ±0.14 
-0.74 ±0.14 


-0.14 ±0.01 
-0.14 ±0.01 
-0.064 ±0.007 
-0.064 ±0.007 


-0.83 ±0.04 
-0.10 ±0.03 
-0.60 ±0.03 
±0.18 ±0.03 



Table 1: Four possible solutions for the model parameters as determined by 
the D+ K*H+iyi data. 
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decay 


B 


[%] 






r+/r_ 


D° ^ K*- 


1.8 ±0.2 
(2.0 ±0.4) 


1.23 ±0.13 


0.16 ±0.04 




1.7 ±0.1 


1.2 ±0.1 


0.16 ±0.04 


(1.88 ±0.29) 


(0.6 ±0.2) 






0.17 ±0.02 


1.34 ±0.2 


0.15 ±0.10 


D+^p" 


0.22 ±0.02 
(< 0.37) 


1.4 ±0.2 


0.15 ±0.10 




0.21 ±0.02 


1.4 ±0.2 


0.16 ±0.10 


Df K*^ 


0.17 ±0.02 


1.3 ±0.2 


0.15 ±0.10 



Table 2: The branching ratios and polarization ratios for the D semilep- 
tonic decays. Where available, the experimental data is quoted in brackets. 



decay 


Bi 


B2 


exp. 




9.4 ±0.5 


9.4 ± 0.5 


6.7 ±0.8 




3±3 


2±2 




Dt v' 


1.6 ±0.7 


0.9 ±0.5 






4±3 


3±3 


7.4 ±3.2 




0.47 ±0.05 


0.5 ±0.5 


0.39^^:^^ 




0.59 ±0.06 


0.7 ±0.6 


0.57 ±0.22 


D+ ^rj 


0.18 ±0.05 


0.1 ±0.2 




D+ rl 


0.021 ±0.005 


0.01 ±0.01 






0.4 ±0.2 


0.2 ± 0.3 





Table 3: The branching ratios for the D ^ P semileptonic decays. 
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decay 


Bth[%] 


^exp 


D+ K*"n+ 


2.4 ± 1.2 


1.92 ±0.19 


D- p-K'^ 


6.6 ± 3.0 


6.6 ± 2.5 


D+ $7r+ 


0.40 ±0.12 


0.61 ±0.06 


L)+ ^ $7r+ 


5.4 ±0.5 


3.6 ±0.9 


Dt P^V 


9.0 ±2.5 


10.3 ±3.2 


Df p+i 


4.5 ± 1.3 


12.0 ±4.5 


D+ X"7r+ 


2.2 ±0.7 


2.74 ±0.29 




4.4 ±0.8 


6.7 ±2.3 




0.029 ±0.005 


0.11 ±0.03 


D+ K*'>p+ 


2.9 ±0.4 


2.1 ± 1.4 


D+ p+r] 


0.05 ± 


< 1.2 


D+ p+r]' 


0.02 ± - 


< 1.5 




0.018 ±0.005 


< 0.28 


D° OUT) 


0.09 ±0.03 




007]' 


0.015 ±0.015 






0.07 ±0.02 


< 0.14 


£)+ ^ $p+ 


0.14 ±0.03 


< 1.5 




0.028 ± 0.004 


< 0.21 



Table 4: The predicted (column two) and measured (column three) branching 
ratios for the nonleptonic decay modes. 
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